BDDC and FETI-DP under Minimalist Assumptions 



Jan MandeH Bedfich Sousedik*™ 

February 1, 2008 



Abstract 

The FETI-DP, BDDC and P-FETI-DP preconditioners are derived in a particulary simple 
abstract form. It is shown that their properties can be obtained from only on a very small set of 
algebraic assumptions. The presentation is purely algebraic and it does not use any particular 
definition of method components, such as substructures and coarse degrees of freedom. It is 
then shown that P-FETI-DP and BDDC are in fact the same. The FETI-DP and the BDDC 
preconditioned operators are of the same algebraic form, and the standard condition number 
bound carries over to arbitrary abstract operators of this form. The equality of eigenvalues of 
BDDC and FETI-DP also holds in the minimalist abstract setting. The abstract framework is 
explained on a standard substructuring example. 



1 Introduction 

The BDDC and FETI-DP methods are iterative substructuring methods that use coarse degrees 
of freedom associated with corners, edges or faces between subdomains, and they are currently the 
most advanced versions of the BDD and FETI families of methods. The BDDC method introduced 
by Dohrmann [3] is a Neumann-Neumann method of Schwarz type [I]. The BDDC method iterates 
on the system of primal variables reduced to the interfaces between the subdomains, and it can be 
understood as further development of the BDD method by Mandel [17]. The FETI-DP method 
by Farhat et al. [5j [6] is a dual method that iterates on a system for Lagrange multipliers that 
enforce continuity on the interfaces, and it is a further development of the FETI method by Farhat 
and Roux [7J. Algebraic relations between FETI and BDD methods were pointed out by Rixen 
et al. [23], Klawonn and Widlund [12j . and Fragakis and Papadrakakis A common bound on 
the condition number of both the FETI and the BDD method in terms of a single inequality in 
was given in [12J. In the case of corner constraints only, a method same as BDDC was suggested 
by Cros [2] . Fragakis and Papadrakakis [9] derived primal versions of FETI and FETI-DP, called 
respectively P-FETI and P-FETI-DP, and they have also observed that the eigenvalues of BDD 
and a certain version of FETI are identical. Mandel, Dohrmann, and Tezaur [19] have proved that 
the eigenvalues of BDDC and FETI-DP are identical and they have obtained a simplifed and fully 
algebraic version (i.e., with no undetermined constants) of a common condition number estimate 
for BDDC and FETI-DP, similar to the estimate by Klawonn and Widlund [12] for BDD and 
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FETI. Simpler proofs of the equality of eigenvalues of BDDC and FETI-DP were obtained by Li 
and Widlund [15], and by Brenner and Sung [I], who also gave an example when BDDC has an 
eigenvalue equal to one but FETI-DP does not. The proof of the equality of eigenvalues of BDD 
and a certain version of FETI and of FETI-DP and P-FETI-DP was recently given by Fragakis [8] 
in a more general framework. 

In this contribution, we derive the FETI-DP, BDDC and P-FETI-DP preconditioners in a 
particulary simple abstract form, with only a very small set of algebraic assumptions (Sec. [3]). 
The presentation is purely algebraical and it does not use any particular definition of method 
components, such as substructures and coarse degrees of freedom. We show that P-FETI-DP and 
BDDC are in fact the same. We then present the condition number bound and the proof of the 
equality of eigenvalues of BDDC and FETI-DP, in the minimalist abstract setting (Sec. E|). Finally, 
we illustrate the abstract framework on a substructuring example (Sec. [5|). 

2 Notation and Preliminaries 

All spaces in this paper are finite dimensional linear spaces. The dual space of a space V is denoted 
by V and (•, •) is the duality pairing. For a linear operator L : W — > V we define its transpose 
L T : V' -» W by (v, Lw) = (L T v,w) for all v G V', w G W, and \\v\\ K = (Kv, v) denotes 
the norm associated with a symmetric and positive definite operator K : V — ► V, i.e., such that 
(Kv, v) > for all v G V, v ^ 0. The norm of a linear operator E : V — ► V subordinate to this 
vector norm is defined by = max„ e v>^o \\Ev \\ K / \\v\\ K . The notation Iy denotes the identity 

operator on the space V. 

Mappings from a space to its dual arise naturally in the variational setting of systems of linear 
algebraic equations. An an example, consider an n x n matrix A and the system of equations 
Ax = b. The variational form of this system is 

x G V : (Ax,y) = (b,y) Vy G V, 

where V = M. n and (•, •) is the usual Euclidean inner product on M. n . For a fixed x, instead of 
the value Ax, we find it convenient to consider the linear mapping y >— > (Ax,y). This mapping is 
an element of the dual space V. Denote this mapping by Kx and its value at y by (Kx,y); then 
K : V — ► V' is a linear operator from V to its dual that corresponds to A. This setting involving dual 
spaces is convenient and compact when dealing with multiple nested spaces, or with dual methods 
(such as FETI). Restricting a linear functional to a subspace is immediate, while the equivalent 
notation without duality requires introducing new operators, namely projections or transposes of 
injections. Also, this setting allows us to make a clear distinction between an approximate solution 
and its residual, which is in the dual space. It is beneficial to have approximate solutions and 
residuals in different spaces, because they need to be treated differently. 
We wish to solve a system of linear algebraic equations 

Ku = f, 

where K : V — > V' , by a preconditioned conjugate gradient method. Here, a preconditioner is a 
mapping M : V — > V. In iteration k the method computes the residual 

r (k) = Ru (k) _f €V ', 

and the preconditioner computes the increment to the approximate solution as a linear 
combination of the preconditioned residual Mr^ k ' G V with preconditioned residuals in earlier 
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iterations. Convergence properties of the method can be established from the eigenvalues A of the 
preconditioned operator MK; the condition number 

= A max (Mif) 
K X min (MK) ' 

gives a well-known bound on the error reduction, cf. e.g. |10j . 

|| e W|| <2 ( - V |lc (0) ll 
where e( fe ) = — u is the error of the solution in iteration k. 

3 Abstract Formulation of the Preconditioners 

3.1 Setting and Assumptions 

We now list a minimalist set of spaces, linear operators, and assumptions needed to formulate the 
BDDC and the FETI-DP methods and to prove their properties. Let TV be a finite dimensional 
space and let a(-,-) be a symmetric positive semi-definite bilinear form on W. Let W C W be a 
subspace such that a is positive definite on W, and / G W. We wish to solve a variational problem 

ueW: a (u,v) = (/,«) Vuef. (1) 

The preconditioners we are interested in are characterized by a selection of an intermediate 
space W, 

W C W C W, (2) 

such that 

a (-, •) is positive definite on W, (3) 
and a selection of linear operators E, B, and Bp. The operator E is a projection onto W, 

E:W^W, E 2 = E, ranged = (4) 
The role of the operator B is to enforce the condition u £ W by Bu = 0, 

5:W^A, nullB = W, ranged = A. (5) 
The operator B^ is a generalised inverse of B, 

B^-.A^W, BB T D = I A . (6) 

The properties ^ - f6]j are enough for the BDDC and the FETI-DP theories separately, and they 
will be assumed from now on. To relate the two methods, we shall also assume that 

B T D B + E = I (7) 

when needed. No further assumptions are made in the rest of the paper. 

To formulate the preconditioners and their properties, we need to define several more linear 
operators from the concepts already introduced. Denote by 

R : W -> W, R : w G W i — >w£W, (8) 

the natural injection from W to W. Clearly, 

ER = I W . (9) 
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Remark 3.1. The conditions (0j - (T^ are satisfied in the applications of FETI-DP and BDDC 
methods, see flStf and Sec. Note that the assumption |2J) allows the case when B is a matrix 
that does not have full row rank. All that is needed is to define A as range -B. In the literature, 
J22J/ and references therein, the projection E is often written in the form E = RDpR T where R 
is a mapping of another space (isomorphic to W) into W . In the abstract setting here, we choose 
to formulate the methods in the space W directly, it turns out that the space W is not needed 
for the theory at all, and R becomes the identity embedding of W into W . The equation {6]) is 
found already in [23, Lemma 1] in a special case. It was extended to form used presently and to 
cover more general algorithms, and used to obtain important connections between dual and primal 
sub structuring methods in \12$ . 

We also define the linear operators S and S associated with the bilinear form a on the spaces 
W and W, respectively, by 

S:W^W', (Sv,w\ =a(v,w) Vv,w €W. (10) 

S:W->W' t (Sv,w^ =a(v,w) \/v, w G W, (11) 

From pop , the variational problem ((TJ) becomes 

Su = f. (12) 
Further, it follows from (TDK, (JIT), and that 

S = R T SR. (13) 

3.2 BDDC 

The following presentation of BDDC follows [20] . It is essentially same as the approach of [1], 
and related to the concept of subassembly in [15J. The BDDC [3\ is the method of preconditioned 
conjugate gradients applied to the system (fT2|) . with the abstract BDDC preconditioner Mbddc '■ 
W — > W defined as 

Mbddc '-r* — >u = Ew, w G W : a(w,z) = (r,Ez) , Vz G W. (14) 

For the equivalence of (j!4[) with other formulations of BDDC, see [19] Lemma 7]. 

From the definitions of S in (|lip and R in (|8|), it follows that the operator form of the BDDC 
preconditioner is 

Mbddc = ES^E?. (15) 

3.3 FETI-DP 

This presentation of FETI-DP follows [19J. The variational problem ([1]) is equivalent to the 
minimization ^ 

—a (u, u) — (/, u) — > min subject to u G W. (16) 

Using 

(f,u} = (f,Eu} = (E T f,u), ueW, 
we can write (I16p as a constrained minimization problem posed on W, 

—a (u, u) — (E T f, u) — > min subject to u G W and Bu = 0. 
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Introducing the Lagrangean 

C(w, A) =~a (u, u) - (E T f, u) + (B T X, u) , 

where A € A' are the Lagrange multipliers, we obtain that problem is equivalent to solving the 
saddle-point problem [22] 

min max£(w, A). (17) 

W (Z\V AeA ' 

Since 

min max£(ii;, A) = max min C(w, A), 
w£ w AeA ' AeA ' wGW 



0, (18) 



it follows that (|16|) is equivalent to the dual problem 

ar(A) _ 
d\ 

where 

^*(A) = mm C(w, A). (19) 
Problem (|18p is equivalent to stationary conditions for the Lagrangean C, 

^£(w,X)lW, (20) 

which is the same as solving for w £ W and A £ A' from the system 

Su, + £ T A = £ T /, , , 

Bw = 0. [ZL) 

Solving from the first equation in (|21|) and substituting into the second equation, we get the dual 
problem in an operator form, 

BS- l B T X = BS^E 7 /. (22) 

The FETI-DP method [6] is the method of preconditioned conjugate gradients applied to the 
problem (|22|) . with the preconditioner given by 

Mfeti-dp = BdSBJ). (23) 

The FETI-DP method solves for the Lagrange multiplier A. The corresponding primal solution 
is found as the minimizer of w in (|19p . or equivalently, from (|20p . which is the same as the first 
equation in (f2"Tj) ; hence, 

w = S- 1 (E T f - B T X) . 

If A is the exact solution of the dual problem (fT8|) , then w GW and so u = w is the desired solution 
of the primal minimization problem (|16p . However, for approximate solution A, in general w ^ W, 
and so the primal solution needs to be projected onto W. We use the operator E for this purpose. 
So, for an arbitrary Lagrange multiplier A, the corresponding approximate solution of the original 
problem is 

u = ES' 1 (E T f - B T X) . (24) 

Note that the operator E does not play any role in FETI-DP iterations themselves. It only serves 
to form the right-hand side of the constrained problem (|21l) . and to recover the primal solution. 
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3.4 P-FETI-DP 



The P-FETI-DP preconditioner [9] is based on the approximate solution from the first step of 
FETI-DP, starting from A = 0, and with the residual r as the right-hand side. The primal solution 
corresponding to the result of this step is the output of the preconditioner. Thus, from (|24p with 
A = and / = r, we have 

M P _ FET i-Dpr = ES~ l E T r. (25) 
Comparing (|25p with the BDDC preconditioner in (|15p . we have immediately: 

Theorem 3.1. The P-FETI-DP and the BDDC preconditioners are the same. 

4 Condition Number Bounds and Eigenvalues 

From (I13h and (|15p . the preconditioned operator of the BDDC method is 

Pbddc = (ES-^EF) (R T SR) . (26) 

From (|22p and (|23|) . the preconditioned operator of FETI-DP is 

Pfeti-dp = (B D SB£) [BS- l B T ^ . (27) 

Clearly, both the BDDC and the FETI-DP preconditioned operators have the same general 
form 

(LA~ l L T ) (T t AT) , (28) 

where A = S is symmetric, positive definite, and L and T are some linear operators such that 

LT = I, (29) 

because of ([6]) and ([9]). This important observation was made in [15] in the equivalent form 
PS^^^PS : range P — > range P, where P is a projection, and in the present form in [1]. 

4.1 Results for the Abstract Form of the Preconditioned Operators 

It is interesting that the fundamental eigenvalue estimate can be proved for arbitrary operators of 
the form (|28p - (|29p . The following lemma was proved in terms of the BDDC preconditioner in [19} 
Theorem 25], and the proof carries over. Because the translation between the two settings is time 
consuming, the proof (with some simplifications but no substantial differences) is included here for 
completeness. The resulting proof of the condition number bound for FETI-DP in Theorem 14.11 
below appears to be new. 

Lemma 4.1. Let V and U be finite dimensional vector spaces and A : V — > V be an SPD operator. 
If L : V — > U and T : U — > V are linear operators such that LT = I on U, then all eigenvalues A 
of the operator (Lyl _1 L T ) (T T AT) satisfy 

1< \ <\\TL\\ 2 A . (30) 
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Proof. The operator (LA 1 L T ) (T T AT) is selfadjoint with respect to the inner product 
(T T ATu,v). So, it is sufficient to bound ((LA~ 1 L T ) (T T AT) u, u) in terms of ( (T T AT) u, u). 
Let u G U. Then 

(LA~~ 1 L T ) (T t AT) u = Lw, (31) 

where w = A~ X L T T T ATu satisfies 

weV, (Aw, v) = (T T ATu, Lv) G V. (32) 
In particular, from (|32p with v = w and (|3ip 

(^w, u>) = (T T ATu, Lw) = (T T ATu, (LA~ l L T ) (T T AT) u) . (33) 
and using LT = /, (f32l) with v = Tu, Cauchy inequality, the definition of transpose, and (f33l) . 



(T T ATu,u) 2 = (T T ATu, LTuf 
= (Aw,Tu) 2 
< (Aw, w) (ATu, Tu) 
= (Aw,w) (T T ATu,u) 

= (T T ATu, (LA~ 1 L T ) (T t AT) u) (T t ATu,u) . 
Dividing by (T T ATu,u) , we get 

(T T ATu,u) < (T T ATu, (LA~ l L T ) (T T AT) u) Vu G U, 



which gives the left inequality in (|30j) . 

To prove the right inequality in (|30p . let again u £ U. Then, from (|3ip . Cauchy inequality in 
the T T AT inner product, definition of the A norm, properties of the norm, and (|33D . 

(T T ATu, (LA~ l L T ) (T T AT) uf 
= (T T ATu,Lw) 2 

< (T T ATu,u) (T T ATLw,Lw) 
= (T T ATu, u) (ATLw, TLw) 
= (T T ATu,u) \\TLw\\ 2 A 

< (T T ATu,u) \\TL\\ 2 A \\wf A 

< (T T ATu,u) \\TL\\\ (T T ATu, (LA~ 1 L T ) (T T AT) u) . 
Dividing by (T T ATu, (LA- 1 L T ) (T t AT) u), we get 

(T T ATu, (LA~ 1 L T ) (T t AT) u) < \\TL\\\ (T T ATu, u) Vn G U. 

□ 

The lower bound in Lemma 14.11 was also proved in a different way in [H Lemma 3.4]. The 
next abstract lemma is the main tool in the comparison of the eigenvalues of the preconditioned 
operators of BDDC and FETI-DP. 
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Lemma 4.2 ([H Lemmas 3.4 - 3.6]). Let V and Ui, i = 1, 2, be finite dimensional vector spaces 
and A : V — > V' be an SPD operator. If Li : V — ► U and Tj : Ui — > V are linear operators such 
that 

LiTi = L onU, i= 1,2, (34) 
T\L\ + T 2 L 2 =1 onV, (35) 

then all eigenvalues (except equal to one) of the operators (LiA~ l Lj^j (T^AT-A and 
{T2 AT2) {L^A^I^) are the same, and their multiplicities are identical. 

4.2 Results for FETI-DP and BDDC 

Condition number bounds now follow immediately from Lemma 14.11 

Theorem 4.1. The eigenvalues of the preconditioned operators of FETI-DP and BDDC satisfy 
1 < A < ojfeti-dp and 1 < A < lobddc, respectively, where 

UBDDC = \\E\\~, UJFETI-DP = \\B^Bw\\ 2 ~. (36) 
In addition, ifW^W and ([?]) holds, then also 

ubddc = ^feti-dp- (37) 

Proof. The eigenvalue bounds with (|36p follows from the form of the preconditioned operators 
()26p and (I27p and Lemma 14.11 The equality (I37p follows from the fact that E and Bj^B are 
complementary projections by ([7|), and the norm of a nontrivial projection depends only on the 
angle between its range and its nullspace [Tl]. □ 

The result in Theorem 14. II was proved in a different way in [18] for BDDC and in [21] for FETI- 
DP. For a simple proof of the bound for BDDC directly from the variational formulation (|14|) . see 
[201 Theorem 2}. 

Equality of the eigenvalues of the two methods follows immediately from Lemma 14.21 

Theorem 4.2. Let f3J) - (0) hold. Then, (a) the spectra of the preconditioned operators of BDDC 
and FETL-DP are the same except possibly for eigenvalue equal to one, and all eigenvalues are 
larger or equal to one, and (b) the multiplicity of any common eigenvalue different from one is the 
same, and the multiplicity of the eigenvalue equal to one for FETI-DP is less than or equal to the 
multiplicity for BDDC. 

Statement (a) of Theorem 14.21 was proved in [19J in a different way, and an elegant simplified 
proof was given in [15] . Statement (b) was proved in [1] . This presentation uses the fundamental 
lemma and the approach from [I]. 

5 Substructuring for a Model Problem 

To clarify ideas, we show how the spaces and operators arise in the standard substructuring theory 
for a model problem obtained by a discretization of the second order scalar elliptic problem. 
Consider a bounded domain $7 C M. d decomposed into nonoverlapping subdomains Qi, i = 1, N, 
which form a conforming triangulation of the domain 0. Each subdomain Oj, from now called 
a substructure, is a union of Lagrangean PI or Ql finite elements, and the nodes of the finite 
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Figure 1: Schematic drawing of continuity conditions between substructures, in the case of corner 
coarse degrees of freedom only: all degrees of freedom continuous (the space W), only the coarse 
degrees of freedom need to be continuous (the space W), and no continuity conditions (the space 
W). 



elements between the substructures coincide. The nodes contained in the intersection of at least 
two substructures are called boundary nodes. The union of all boundary nodes of all substructures 
is called the interface, denoted by T, and Tj is the interface of substructure Qj. The space of 
all vectors of local degrees of freedom on Tj is denoted by Wj. Let Si : Wi — > Wj be the Schur 
complement operator obtained from the stiffness matrix of the substructure fij by eliminating all 
interior degrees of freedom of Oj, i.e., those that do not belong to interface IV We assume that the 
matrices S% are symmetric positive semidefinite. Let W = W\ x • • • x Wn and write vectors and 
matrices in the block form 



w 



Wl 



W N 



Si 



The bilinear form a is then given by 



w 6 W, S 



a (u, v) = u T Sv. 



Sn 



(38) 



The solution space W of the problem ([I]) is a subspace of W such that all subdomain vectors of 
degrees of freedom are continuous across the interfaces, which here means that their values on all 
the substructures sharing an interface nodes coincide. 

The BDDC and FETI-DP are characterized by selection of coarse degrees of freedom, such as 
values at the corners and averages over edges or faces of subdomains (for their general definition see, 
e.g., |13j). In the present setting, this becomes the selection of the subspace W C W, defined as the 
subspace of all functions such that coarse degrees of freedom are continuous across the interfaces. 
Cf., Fig. [TJ There needs to be enough coarse degrees of freedom that the variational problem on W 
is coercive, i.e., (|3|) is satisfied. Creating the stiffness matrix on the space W is called subassembly 

The last ingredients are the selections of the linear operators E, B, and Bq. The operators E 
and B are in fact defined on the whole space W; they are considered restricted on W only for the 
purposes of the theory here. The operator E : W — * W is an averaging of the values of degrees 
of freedom between the substructure. The averaging weights are often taken proportional to the 
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diagonal entry of the stiffness matrices in the substructures. The matrix B enforces the continuity 
across substructure interfaces by the condition Bw = 0. Each row B has only two nonzero entries, 
one equal to +1 and one equal to —1, corresponding to the two degrees of freedom whose value 
should be same. So, Bw is the jump of the value of w between substructures. Redundant Lagrange 
multipliers are possible; then B does not have full row rank and A = range B is not the whole 
euclidean space. Finally, Bq is a matrix such that a vector A of jumps between the substructures 
is made into a vector of degrees of freedom Bj^X that exhibits exactly those jumps. That is, 
BB p = I. The construction of Bp involves weights, related to those in the operator E, so that 
BpB + E = I, and its details are outside of the scope of this paper. Such construction was 
done first for the FETI method in [14] in order to obtain estimates independent of the jump of 
coefficients between substructures, and then adopted for FETI-DP. We only note that in many 
cases of practical relevance, the matrix Bd is determined from the properties (@|) - ([7]) uniquely as 
the Moore-Penrose pseudoinverse in a special inner product given by the averaging weights in the 
operator E |19[ Theorem 14]. 
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